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OPTIMIZATION  OF  STEEL  PILE  FOUNDATIONS  USING  OPTIMALITY 


CRITERIA 


1 .  INTRODUCTION 

Optimization  techniques  have  come  under  greater  use  as 
the  demand  has  grown  for  least  cost  structures.  Modern 
structural  optimization  began  to  be  developed  in  the  1960's 
(1,2,3).  Two  categories  of  optimization  methods  exist: 
analytic  and  numerical  methods  (4,5).  Analytic  methods  can 
lead  to  the  exact  optimal  solution,  but  these  methods  require 
all  of  the  optimization  parameters  to  be  expressed  explicitly 
in  terms  of  the  optimization  variables.  Numerical  methods  are 
used  for  structural  optimization  problems.  They  are  used 
because  the  problem  does  not  need  to  be  explicitly  expressed. 
Numerical  methods  usually  involve  iterative  processes  which 
alter  an  initial  design  towards  an  optimal  design.  Each 
optimization  method  has  advantages  and  disadvantages.  The 
optimization  method  must  be  selected  after  considering  the 
characteristics  of  the  problem  to  be  solved. 

The  layout  of  piles  is  to  be  optimized.  The  objective  of 
the  optimization  is  to  find  the  least  weight  of  steel  HP-14 
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to  the  design  variables  which  can  be  written  as: 


dd,  fy  dd t 


=  0 


i=l,  .  .  .  ,n 


(2.8) 


where  there  are  n  design  variables  d;.  This  is  the  first 
necessary  condition.  Note  that  Xj  does  not  vary  at  the 
optimum. 

The  other  necessary  conditions  are  found  by  looking  at 
active  and  passive  constraints.  A  constraint  gradient  is  used 
in  the  vector  sum  to  oppose  improvements  in  weight  only  when 
the  constraint  is  active.  This  requires  that  Xj  is  zero  when 
hj  is  not  zero,  and  Xj  is  not  equal  zero  when  hj  is  zero. 

kjhj  -  0  j- 1,  .  .  .  ,m  (2.9) 

The  final  conditions  dictate  that  the  Lagrange  Multipliers 
must  be  positive  and  the  constraint  must  not  be  violated. 

^  i  0  j»l, . . . ,m  (2.10) 


hj  s  0  j*  1,  .  .  .  ,m  (2.11) 

Selecting  the  correct  set  of  constraints  as  active  is  the  most 
complicated  portion  of  the  optimization  process.  Equation  2.8 
may  be  rewritten  to  provide  the  optimality  criteria: 


(2.12) 


This  criteria  must  be  satisfied  for  the  n  variables. 


2 


The  optimality  criteria  for  the  i*  variable  is  used  as 
the  efficiency  of  the  variable.  The  i*  variable  (d;)  is 
altered  by  the  value  of  the  i*  optimality  criteria  in  the 
following  recurrence  formula: 


where  v  is  the  index  of  the  iteration  number  and  r  is  the 
convergence  control  parameter. 

The  optimality  criteria  is  less  than  one  when  the 
component  of  the  objective  or  weight  gradient  for  the  i* 
variable  is  greater  than  the  resistance  provided.  When  the 
weight  gradient  is  positive,  the  variable  is  free  to  decrease 
by  an  amount  roughly  proportional  to  the  value  of  the 
optimality  criteria.  If  the  resistance  is  stronger  than  the 
weight  gradient  then  the  optimality  criteria  is  greater  than 
one  and  the  variable  should  be  increased.  Taking  the  r"1  root 
of  the  optimality  criteria  where  r>l  assures  that  the 
prediction  of  d;  for  the  next  iteration  does  not  greatly 
overshoot  the  optimum  value.  A  reasonable  value  for  the 
convergence  control  parameter  (r)  is  taken  as  2. 

Rewriting  Equation  2.13: 


(2.14) 
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and  expanding  it  using  the  binomial  theorem  results  in  a 
linear  recurrence  eguation: 


d/*1  =  dr 


i+— 

r 


'y  x.^i) /*) 

\M  ’ddjddj 


(2.15) 


The  values  of  the  Lagrange  Multipliers  to  be  used  in  the 
variable  recurrence  relation  must  be  estimated.  Estimation  of 
the  multipliers  may  be  performed  using  either  a  recurrence 
relationship  or  by  solving  a  set  of  linear  equations.  A 
recurrence  relation  method  requires  assuming  the  initial 
values  of  the  Lagrange  Multipliers.  Another  disadvantage  is 
that  the  process  is  slow  to  converge.  The  second  method  which 
was  the  method  used  was  a  linear  equation  solution  method. 
This  method  requires  that  the  constraints  must  be  known  a 
priori  to  be  contained  in  either  the  active  or  passive  set  of 
constraints. 

A  formula  for  the  change  in  the  active  constraints  is 
used  to  produce  the  linear  equation  solution  method.  The 
equation  for  the  change  and  linear  approximations  of  the 
change  in  the  j*  active  constraint  can  be  written  as: 

n  gu  (2.16) 

Afy-  hjidt+bdt)  ~hj(dj)  =  ^  ~dd.^di 


where  hj(dj)  is  the  value  of  the  constraint  given  the  variable 
dj.  The  expected  value  of  the  active  constraints  after  any 
iteration  should  approach  zero,  therefore  the  previous  value 
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sections  that  can  be  used  to  satisfy  the  demands  of  the 
loading  conditions. 

Previous  work  has  been  performed  in  the  area  of  pile 
optimization.  A  program  was  developed  in  1981  by  Hill  which 
attempts  to  optimize  the  layouts  of  piles  (6) .  To  reduce  the 
weight  of  steel  a  pile  deletion  process  was  used.  The 
optimization  procedure  consisted  of  first  finding  the  optimal 
pile  slopes,  then  finding  the  optimal  pile  spacing  within 
specified  zones,  and  deleting  piles  until  the  stresses  and 
displacements  are  near  their  limits.  The  piles  are  deleted  in 
an  iterative  process  by  eliminating  the  most  and/or  the  least 
stressed  piles. 

An  optimality  criteria  approach  was  used  in  this  study. 
The  pile  optimization  problem  is  nonlinear.  The  design  limits 
such  as  stresses  change  nonlinearly  with  respect  to  the 
variables.  Frequently  problems  are  linearized  by  finding  the 
gradients  with  respect  to  the  variables  of  the  weight  and  the 
constraints. 

The  optimality  criteria  approach  is  developed  from  the 
nature  of  an  optimal  design  point.  Optimality  criteria 
methods  were  first  used  in  the  1970's  and  have  been  applied  to 
a  variety  of  structural  optimization  problems.  A  paper  by 
Venkayya  presents  the  optimality  criteria  method  based  on  a 
strain  energy  criterion  (7)  .  The  method  is  developed  to 
solve  optimization  problems  with  a  large  number  of  variables. 
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The  method  quickly  arrives  at  an  optimal  design  while  using  a 
relatively  small  amount  of  computational  effort.  Venkayya 
applies  the  method  to  truss  and  frame  optimization  by  allowing 
member  sizes  to  vary. 

The  optimization  methods  which  are  developed  are  applied 
to  an  existing  pile  group  analysis  computer  program.  See 
Appendix  1  for  further  information  about  the  analysis  methods. 
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2.  MATHEMATICAL  FORMULATION  OF  THE  METHOD 

2.1  STATEMENT  OF  THE  FORM  OF  THE  OBJECTIVE  AND  CONSTRAINTS 

The  objective  function  of  the  optimization  process  is  the 
function  to  be  minimized.  The  goal  of  the  pile  optimization 
process  is  to  find  the  minimum  weight  of  steel  that  can  be 
used  to  satisfy  the  given  constraints.  The  objective  function 
is: 

wTM£Pivi  <2-*> 

where  Wx  is  the  total  weight,  P}  is  the  density,  and  V;  is  the 

volume  for  the  structural  element  i,  and  N  is  the  number  of 
piles. 

Structures  are  subject  to  a  variety  of  constraints.  Some 
of  the  constraints  are:  member  stress  limits,  displacement 
limits,  soil  capacity,  minimum  and  maximum  member  sizes,  and 
the  layout  restrictions.  The  constraints  generally  place  an 
upper  or  lower  bound  on  the  value  of  a  parameter  such  as 
stress.  An  example  of  a  stress  constraint  is: 

hj  =  Oj-'Oji  0  j- 1,...,%  (2.2) 

where  hj  is  the  value  of  the  constraint  j,  a}  is  the  stress, 

is  the  upper  bound  of  the  stress  in  element  j,  and  m,  is  the 
number  of  stress  constraints.  The  constraint  is  satisfied 
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when  hj  is  zero  or  negative.  To  improve  the  numerical 
conditioning  of  the  constraints,  Equation  2.2  may  be  rewritten 
as: 

o  • 

hj  -  -=^  -  1  SO  J=l,  .  .  .  ,11^  (2.3) 

A  constraint  on  the  lower  bound  of  a  parameter  such  as  member 
size  can  be  written  as: 

hj  '=  ZlzLl  +  i  s  o  j*l ,  .  .  . , (2.4) 

j 

where  X„j  is  the  lower  bound  of  the  moment  of  inertia  IMj  for 
member  j,  and  m2  is  the  number  of  lower  limit  size 
constraints. 

2.2  OPTIMALITY  CRITERIA  FORMULATION 
The  problem  statement  is: 
minimize: 

<*-S> 

subject  to: 

hj  s.  0  j= l,  .  .  .m  (2.6) 

where  m  is  the  total  number  of  constraints. 

Satisfaction  of  the  KUHN-TUCKER  conditions  is  necessary 
for  a  point  in  design  space  to  be  a  local  minimum  of  the  given 
problem  (8)  .  One  optimal  point  exists  when  a  problem  is 
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convex.  If  the  design  space  is  convex  then  the  satisfaction 
of  the  KUHN-TUCKER  conditions  is  sufficient  for  a  design  to  be 
the  global  optimum.  The  pile  orientation  optimization  problem 
is  non-convex.  A  method  that  can  be  used  to  find  the  global 
optimum  is  the  usage  of  multiple  starting  points  for  the 
optimization  process.  Forming  the  KUHN-TUCKER  conditions 
requires  the  Lagrangian  which  is: 

L  *  wt  +  £  * jhj  (2.7) 
j*  i 

where  Xj  is  the  Lagrange  Multiplier  for  the  constraint  j. 

The  first  necessary  condition  for  a  local  minimum 
dictates  that  a  further  improvement  in  the  weight  is  prevented 
by  constraints.  See  Figure  2-1.  Constraints  are  active  when 
h—0,  and  passive  if  hj<0.  The  constraints  which  are  active  at 
the  optimum  design  point  will  form  a  vector  opposing  the 
weight  gradient.  In  Figure  2-1  constraints  1  and  2  are 
active.  To  perform  a  vector  addition  of  weight  and  constraint 
gradients  a  scaling  variable  is  required  to  account  for  the 
units  of  the  gradients.  The  scaling  is  performed  by  the 
Lagrange  Multiplier.  The  formula  for  the  vector  addition  of 
the  gradients  is  the  derivative  of  the  Lagrangian  with  respect 
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(a) 


Weight  vector 


Gradient  of  constraint  1.  Gradient  of  constraint  2. 

(b) 


Figure  2-1  Gradient  Vectors 

(a):  Constrained  space,  (b):  Gradients. 
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of  the  constraint  should  be  equal  to  the  negative  of  the 
change . 


-  - 1 Is2 


j  =  l*  • 


/in 


(2.17) 


The  change  in  the  design  variables  can  be  found  using  the 
variable  recurrence  equation  which  contains  the  first  KUHN- 
TUCKER  necessary  condition.  Substituting  the  variable  change: 


A  =  +1-  (2.18) 

and  Equation  2.15  into  the  constraint  change  equation  yields: 


7  f  ,  dh  BWT)  " 
V  k  ***/  Mi) 


(2.19) 


d/  j  =1, . . ./n 


and  then  rearranging  yields: 


^  dhs  dha  ,  dWT\^  A  3A,  '*• 

hj  (di)  S  AsS(^  3dj 1  adjdi+  S  7  1 . 777 


(2.20) 


Passive  variables  are  variables  which  are  forced  to 
maintain  a  certain  value.  A  cross-section  size  can  be  held  at 
the  maximum  or  minimum  value  by  making  it  passive.  An 
additional  term  is  required  to  account  for  the  change  in  the 
constraints  by  the  applied  change  in  the  passive  variables. 
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....  A.£ Jdh,  dhe  .  dWT\^  Xdh,^  °  dh,  .  . 

^  -SS5N**'  *gp*  £  ^d<* 


j -1, . . . ,m 


(2.21) 


where  n,  is  the  number  of  active  variables. 

2.3  SELECTING  THE  ACTIVE  CONSTRAINT  SET 

The  solution  method  consists  of  assuming  an  active  set  of 
constraints  and  forming  the  linear  equations.  If  the  solution 
of  the  equations  satisfies  the  KUHN-TUCKER  conditions 
(Equations  2.8-2.11)  then  the  solution  is  the  optimum  for  the 
current  iteration.  The  variable  recurrence  equations  can  then 
be  solved. 

The  process  of  selection  of  the  active  Lagrange 
Multipliers  can  be  performed  by  first  choosing  the  constraints 
which  are  near  violation  as  active.  The  set  of  linear 
equations  is  formed  for  the  active  constraints.  The  set  is 
altered  by  causing  newly  violated  constraints  to  be  considered 
active  and  removing  constraints  with  negative  Lagrange 
Multipliers.  To  aide  in  approaching  a  beneficial  set  it  may 
be  decided  that  at  any  trial  only  the  most  violated 
constraints  will  enter  the  solution.  When  beginning  a  new 
optimization  iteration  a  good  initial  selection  is  the  final 
set  of  active  constraints  from  the  previous  optimization 
iteration. 
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If  a  feasible  set  of  active  constraints  can  not  be  found 
then  an  extension  of  the  search  method  is  used.  It  is 
possible  to  grow  a  set.  First  only  one  constraint  is  given 
consideration.  It  is  considered  the  only  constraint  applied 
to  the  design.  After  a  set  of  active  constraints  is  found 
then  the  number  of  constraints  applied  to  the  design  is 
incremented.  The  new  set  of  active  constraints  is  found 
beginning  with  the  previous  set.  This  process  is  continued 
until  all  of  the  constraints  given  are  applied  to  the  design. 

The  process  of  choosing  active  constraints  could  use 
refinement.  The  selection  process  could  be  solved  more 
directly  by  solving  a  second  order  problem.  The  linear 
equations  have  the  first  KUHN-TUCKER  condition  embedded  in 
them.  The  equations  to  be  solved  would  be  the  linear 
equations  with  the  remaining  KUHN-TUCKER  conditions.  This  is 
a  second  order  problem  since  either  h  or  X  must  be  zero.  A 
numerical  method  for  solving  the  set  of  second  order  equations 
relies  on  linearization.  The  linearization  results  in  a  very 
poorly  conditioned  problem. 
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3.  APPLICATION  OF  OPTIMIZATION  METHOD 

The  pile  layouts  were  analyzed  by  an  altered  version  of 
the  US  Army  Corps  of  Engineers  program:  Pile  Group  Analysis 
(CPGA)  Computer  Program  (9) .  The  optimization  process 
requires  the  constraint  gradients  which  are  the  changes  in  the 
constraints  due  to  a  variable  change.  The  gradients  are  found 
by  repetitively  reanalyzing  the  layout.  The  allowable  value 
of  the  stress  constraints  will  change  with  the  design 
variables.  The  allowable  stresses  are  updated  in  accordance 
with  the  US  Army  Corps  of  Engineers  document:  BASIC  PILE  GROUP 
BEHAVIOR  (10) .  See  Appendix  l  for  further  details  on  the 
analysis  of  pile  groups. 

3.1  PILE  ORIENTATION  PROBLEM 

The  layout  of  piles  under  dam  structures,  retaining 
walls,  and  control  towers  is  to  be  optimized.  The  cross- 
sections  and  the  orientations  of  the  piles  will  be  varied. 
The  initial  coordinates  of  the  head  of  the  piles  will  be 
fixed.  The  batter  is  the  orientation  that  will  vary.  The 
batter  of  the  piles  is  the  ratio  of  the  depth  of  the  pile  to 
the  horizontal  distance  between  the  pile  head  and  tip  as  shown 
in  Figure  3-1. 

Variables  in  the  optimization  process  can  be  varied 
either  discretely  or  continuously.  In  a  continuous  process 
the  variables  may  acquire  any  value  within  the  design  space. 
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Theta  ■  0.  degrees  Theta  ■  90.  degrees 

Lower  case  x  and  y  are  the  local  axes 

(c) 


Figure  3-1.  Pile  coordinate  system, 
(a):  Batter,  (b):  Phi,  (c):  Theta. 
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In  a  discrete  process  the  variables  are  allowed  to  acquire 
only  values  contained  within  a  specified  set  of  values. 

The  cross-sections  will  not  be  restricted  to  discrete 
sizes  but  will  be  governed  by  a  continuous  set  of  equations 
which  describe  the  HP-14  sections.  It  is  required  that  one 
variable  should  entirely  represent  the  cross-section  in  the 
optimization  process.  All  the  cross-sectional  properties  will 
be  expressed  in  terms  of  one  primary  variable  which  is  the 
major  axis  moment  of  inertia  (I„) .  The  secondary  variables 
are  the  minor  axis  moment  of  inertia  (Iyy)  ,  the  area  (A)  ,  and 
the  extreme  fiber  distances  (cxf  cy)  .  The  following  equations 
are  approximations  developed  from  the  AISC  Manual  for  HP-14 


sections  (11).  See  Figures  3-2  to  3-5. 

Iyy  =  0.  70  II xx  -  .2200  (inche  4)  (3.1) 

=  .  1  I  •  1  (3.2) 

cx  =  1 . 2110  *  10"3  Ixx  +  12.71  2  ( inche  )  (3.3) 

Cy  —  .  110  *10'4  Ixx  +  14  •  1  (inche)  (3.4) 


The  pile  orientation  problem  is  subject  to  several  types 
of  constraints  such  as  maximum  member  stress  limits, 
displacement  limits,  soil  capacity,  maximum  and  minimum  member 
sizes,  and  the  layout  restrictions.  The  first  layout 
restriction  demands  that  the  tips  of  two  piles  may  not  come 
into  contact  with  each  other.  In  two-dimensions  the  tip 
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Major  Axis  inertia  (Ixx  in4) 


Figure  3-2.  Minor  Axis  Inertia  vs.  Major  Axis  Inertia 
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A  =  0.049153  Ixx 
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500  1000  1500 

Major  axis  inertia  (Ixx  in4) 

Figure  3-3.  Area  vs.  Major  Axis  Inertia 
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Cy  =  6.11  *  10-4  Ixx  +  14.139  in 


Figure  3-4.  Extreme  fiber  distance  in  y— direction  vs.  major  axis  inertia 


Figure  3-5.  Extreme  fiber  distance  in  x-direction  vs.  Major  axis  inertio 


constraint  is  applied  by  constraining  the  coordinate  of  one 
tip  to  be  always  greater  than  the  coordinate  of  the  second 
tip.  Another  restriction  is  that  piles  may  not  be  battered  at 
angles  that  are  far  from  the  vertical  position.  This  is 
called  a  minimum  batter  constraint. 
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3.2  APPLICATION  OF  OPTIMALITY  CRITERIA 

The  optimization  algorithm  needs  the  gradients  of  the 
constraints  and  the  objective  function.  The  gradients  are 
approximated  using  a  finite  difference  method.  Each  of  the 
design  variables  is  altered  to  find  the  resulting  change  in 
the  constraints. 

3.2.1  ENERGY  OBJECTIVE  FUNCTION 

The  stored  energy  resulting  from  global  forces  and 
displacements  will  be  used  as  an  alternate  objective  function 
for  cases  when  the  weight  objective  is  inadequate.  It  is 
stated  as: 

"r  -  piDi  <3-s> 

where  P;  is  the  global  force,  D*  is  the  global  displacement  for 
component  i,  and  WT  is  the  new  artificial  objective  function. 

Piles  that  are  vertical  have  a  weight  gradient  which  is 
zero  for  a  change  in  batter.  A  non-zero  weight  gradient  is 
needed  by  the  algorithm  to  form  the  linear  Lagrange  Multiplier 
equations.  An  arbitrary  objective  which  is  small  relative  to 
the  other  gradient  components  could  be  used.  The  small  weight 
gradient  would  reflect  that  the  variable  has  no  direct  effect 
on  the  weight.  This  arbitrary  gradient  should  be  large  enough 
so  that  the  solution  is  not  numerically  ill-conditioned.  When 
energy  is  taken  as  an  alternate  objective  of  the  vertical 
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piles  then  the  energy  gradients  can  be  scaled  to  a  level  so 
that  they  are  slightly  significant.  The  energy  gradients  of 
the  batters  are  scaled  to  one  percent  of  the  value  of  the 
weight  gradient  for  a  batter  that  is  not  vertical. 

Using  the  energy  as  an  objective  function  for  the 
vertical  piles  provides  a  greater  meaning  than  using  a 
completely  arbitrary  objective  function.  An  energy  objective 
function  is  effectively  a  displacement  objective  function 
since  the  global  forces  are  constant  and  only  displacements 
will  vary.  One  reason  the  energy  objective  has  meaning  is 
because  small  global  cap  displacements  are  desirable.  Another 
reason  that  this  objective  is  proper  is  that  piles  have  less 
displacement  when  they  are  nearly  aligned  with  the  load.  The 
resulting  alignment  would  help  the  pile  to  resist  the  loads 
more  efficiently  by  axial  deformation  rather  than  by  flexure. 

In  problems  where  piles  are  allowed  to  rotate  in  all 
three  dimensions  the  variable  of  rotation  about  the  vertical 
axis  would  not  have  a  weight  change  associated  with  it.  An 
alternate  objective  such  as  an  energy  objective  could  be  used 
for  these  variables.  This  variable  is  not  considered  in  the 
examples . 

3.2.2  TOPOLOGICAL  VARIABLES 

The  variables  which  represent  the  batter  of  the  piles  may 
be  measured  by  using  either  the  angles  of  rotation  or  the 
coordinates  of  the  tip  of  the  piles.  Using  the  tip 


21 


coordinates  as  the  variables  would  simplify  the  pile  tip 
interference  constraints. 

Each  batter  is  given  it's  own  local  coordinate  system. 
In  this  report  the  distance  between  the  local  coordinate 
system  and  the  current  tip  position  will  be  referred  to  as 
slack. 

The  choice  of  the  origin  of  the  coordinate  system  will 
have  an  effect  on  the  convergence  of  each  variable.  Consider 
this  example:  A  variable  which  has  nothing  restricting  it's 
change  will  prefer  to  passively  decrease  itself  to  zero.  If 
the  convergence  control  parameter  (r)  is  equal  to  two  then  the 
variable  would  halve  itself.  If  this  variable  is  a 
topological  variable  with  a  slack  of  80  feet  then  the  tip 
would  move  40  feet,  however  the  pile  behavior  may  be 
predictable  for  a  change  of  only  30  feet.  The  slack  should  be 
chosen  given  the  nature  of  the  variable  and  it's  gradients. 
If  the  constraints  vary  linearly  with  the  variables  then  the 
variables  have  well  behaved  gradients.  This  greater 
predictability  would  allow  an  increased  slack. 

Certain  topological  variables  improve  the  objective 
function  as  their  value  increases.  They  have  a  negative 
weight  gradient.  An  unrestricted  variable  with  a  negative 
weight  gradient  should  increase.  Given  the  unaltered 
algorithm  all  variables  which  are  not  actively  controlled  by 
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a  constraint  would  decrease  towards  the  origin  of  the  local 
coordinate  system.  When  r  is  two  they  would  halve  themselves. 

The  variable  recurrence  formula  appears  to  need 
reformulation.  It  was  developed  using  the  optimality  criteria 
as  a  measure  of  the  efficiency.  See  Equations  2.12  and  2.13. 
When  the  weight  gradient  is  larger  than  the  restraining 
components  the  variable  is  unrestricted;  which  is  true  for  any 
weight  gradient  value.  The  negative  coefficients  of  the 
negative  weight  gradient  cause  the  resulting  optimality 
criteria  value  to  be  greater  than  1  instead  of  less  than  1. 

Two  methods  could  be  used  to  adapt  the  optimization 
method  without  altering  the  recurrence  equations.  The  first 
is  to  simply  reverse  the  sense  of  the  local  coordinate  system 
so  that  the  weight  gradients  are  always  positive.  The  method 
which  is  used  consists  of  applying  negative  slack  while 
maintaining  the  sense  of  the  coordinates.  This  method  is 
chosen  so  that  all  of  the  gradients  do  not  have  to  be  switched 
in  sign  and  a  small  amount  of  computer  time  is  saved.  For 
this  method  when  a  variable  is  not  controlled  and  when  r  is 
two  it  would  halve  itself.  The  negative  variable  would 
increase  towards  the  origin. 
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4 .  NUMERICAL  RESULTS 

The  examples  optimized  the  pile  layouts  for:  one 
retaining  wall,  four  dams,  and  one  control  tower. 

Each  pile  does  not  have  an  independent  cross-section  and 
orientation.  The  piles  are  linked  into  groups  which  share  the 
same  cross-section  and  batter.  This  will  ease  installation  of 
the  piles,  and  it  will  reduce  the  size  of  the  optimization 
problem  to  be  solved. 

TYPES  OF  CONSTRAINTS 

1.  Member  stress  constraints.  The  most  highly  stressed 
member  in  any  group  is  prevented  from  becoming  over-stressed. 

2.  Size  constraints.  The  cross-sections  are  prevented 
from  advancing  beyond  the  maximum  or  minimum  section  sizes 
available. 

3.  Interference  constraints.  Piles  can  not  pass  through 
one  another.  In  two-dimensional  problems  and  the  less 
complicated  three-dimensional  problems  this  constraint  can  be 
established  by  preventing  the  pile  tips  from  meeting. 

4.  Minimum  batter  constraints.  The  piles  may  not  be 
battered  parallel  to  the  ground  surface  or  beyond  a  specified 
angle. 

5.  Displacement  constraints.  Limits  on  displacements  are 
thought  to  aide  problem  convergence  by  preventing  bizarre 
estimates  of  the  optimum  point. 
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In  these  examples  the  most  frequently  active  constraints 
at  the  optimum  are  the  member  stress,  and  member  size 
constraints. 

4.1  EXAMPLE  1. 

Example  1 .  This  example  is  a  retaining  wall.  The 
example's  symmetry  of  geometry  and  loading  allows  a  segment  of 
the  wall  to  be  modeled  as  a  two-dimensional  problem.  The 
segment  has  five  piles  arranged  into  two  groups  as  shown  in 
Tables  4-1,2  and  Figure  4-1. 

The  final  layout  is  found  in  Table  4-3,  and  Figure  4-2. 
This  example  reaches  a  convergent  solution  in  five  iterations 
of  the  optimality  criteria  algorithm.  The  weight  improves 
from  15294.  to  12766.  lbs  as  shown  in  Figure  4-3.  In  the 
first  iteration  the  weight  increases  because  the  original 
layout  has  a  44%  stress  constraint  violation  in  the  second 
load  case. 

The  active  constraints  at  the  optimum  are  the  minimum 
inertias  of  both  groups,  and  the  stresses  in  members  3  and  1 
due  to  load  cases  1  and  2  respectively.  See  Figure  4-4  for 
the  convergence  of  the  stress  constraint  for  member  1. 
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Table  4-3.  Final  layouts  and  sizes  for  example  1. 


Pile 

Group 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

1 

729. 

4.69 

180. 

0. 

2 

729. 

1.04 

0. 

0. 
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Pile  group  1.  Pile  group  2. 

depth  ■  30' 

<b) 


Figure  4-1.  Initial  pile  layouts  for  examples  1., 11,1.2, 
and  1.3.  (a):Plan  view,  (b):Elevation. 
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L  Bit  tit' 
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Y 

(a) 


Arrows  Indloats 
battor  direction 


Group: 

□  1 


Pile  group  1.  Pile  group  2. 

depth  ■  30’ 

(b) 

Figure  4-2.  Final  pile  layouts  for  example  1. 
(a)-.PIan  view,  (b):Elevation. 
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Example  1.1:  This  example  has  the  geometry  of  example  1 , 
but  the  loading  is  much  less.  The  loading  was  arbitrarily 
reduced  as  shown  in  Table  4-4.  \farious  versions  of  example  1 
were  optimized  to  explore  the  behavior  of  the  optimization 
process.  The  loading  for  this  example  was  reduced  to  be  able 
to  view  the  optimization  performance  for  a  design  which  is 
originally  too  conservatively  designed.  An  optimization 
algorithm  should  significantly  reduce  the  weight  for  a 
initially  conservative  design. 

The  final  layout  is  found  in  Table  4-5,  and  Figure  4-5. 
The  weight  improves  from  15294.  to  10924.  lbs  after  15 
iterations  as  shown  in  Figure  4-3.  The  minimum  size  for  the 
two  groups  are  the  only  constraints  that  are  active.  The 
loading  is  so  small  that  the  pile  batters  become  vertical  to 
minimize  the  weight.  The  pile  batters  are  not  required  to  be 
outwardly  positioned  to  be  aligned  with  the  load  vector. 


3Q 
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1 


2 


3 


Arrow*  Indicate 
battar  dlraotlon 


Pile  group  1.  Pile  group  2. 

depth  ■  30* 

(b) 


Figure  4-5.  Final  pile  layouts  for  example  1.1 
(a):Plan  view,  (b):Elevation. 


Example  1.2;  This  example  is  example  1.1  with  a  pile  tip 
interference  limit  of  six  feet.  Interference  constraints  are 
not  frequently  active  in  the  following  examples.  The  minimum 
distance  of  six  feet  between  the  pile  groups  was  expected  to 
force  this  constraint  to  become  active. 

The  final  layout  is  found  in  Table  4-6,  and  Figure  4-6. 
After  12  iterations  a  final  weight  of  10939.  lbs  is  reached  as 
shown  in  Figure  4-3.  Note  the  weight  is  very  close  to  the 
previous  example.  It  is  expected  that  making  constraints  more 
strict  would  increase  the  weight.  The  interference 
constraints  simply  force  the  batters  from  the  vertical 
positions.  The  insignificant  weight  gain  demonstrates 
insensitivity  of  the  weight  to  the  batters. 

The  active  constraints  are  the  same  as  example  1.1  with 
the  addition  of  the  interference  constraint.  The  convergence 
of  the  interference  constraint  is  found  in  Figure  4-7. 


Table  4-6.  Final  layouts  and  sizes  for  example  1.2. 


■gag 

m 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

i 

729. 

20.3 

180. 

0. 

2 

729. 

18.3 

0. 

0. 
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e* 


Pile  group  1.  Pile  group  2. 

depth  ■  30' 

(b) 


Figure  4-6.  Final  pile  layouts  for  example  1.2 
(a).PIan  view,  (b):Elevation. 
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Figure  4-7.  Example  1.2.  Separation  between  piles  3  and  4. 
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Example  1.3;  In  this  example  the  displacements  are 
restricted  to  view  this  type  of  constraint  when  active.  The 
displacement  in  the  x-direction  is  restricted  to  0.10  inches. 
See  the  initial  layout  in  Tables  4-2,4  and  Figure  4-1. 

The  final  layout  is  found  in  Table  4-7,  and  Figure  4-8. 
After  14  iterations  the  weight  converges  to  11211.  lbs.  The 
weight  increased  slightly  compared  to  example  1.1  as  shown  in 
Figure  4-3.  The  displacement  constraint  in  addition  to  the 
constraints  of  example  l.i  are  active.  The  displacement  is 
precisely  controlled.  The  displacement  convergence  is  found 
in  Figure  4-9. 


Table  4-7.  Final  layouts  and  sizes  for  example  1.3. 


1  Pile 

|  Group 

(in1) 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

1  1 

729. 

13.7 

180. 

0. 

.iHBMD 

729. 

2.81 

0. 

0. 
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Pile  group  1.  P*l®  9rouP  2. 

depth  •  30' 

(b) 


Figure  4-8.  Final  pile  layouts  for  example  1.3 
(a):Plan  view,  (bhElevation. 


Displacement 
Component  1 


01  2  3  4  5  6  7  8  9  10  11  12  13  14 

Iteration 

Figure  4-9.  Satisfaction  of  displacement  constraint  for  example  1.3. 


Example  1.4:  The  choice  of  linking  piles  into  groups  is 
performed  by  the  designer.  It  should  be  determined  if  the 
best  grouping  has  been  chosen.  In  this  version  of  example  1 
the  piles  are  unlinked.  Their  final  layout  should  give 
insight  into  the  best  grouping.  The  pile  tips  are  constrained 
to  be  separated  by  at  least  two  feet. 

See  the  initial  and  final  designs  in  Tables  4-1, 8, 9  and 
Figures  4-10,11.  After  11  iterations  the  weight  has  converged 
from  12799.  to  11753.  lbs  as  shown  in  Figure  4-12.  As 
expected,  this  is  slightly  smaller  than  the  12766.  lbs.  for 
example  1.  with  linked  piles.  The  active  constraints  are: 
the  minimum  size  for  all  groups;  stress  in  piles  1  and  2  for 
load  cases  1  and  2  respectively;  and  the  interferences  of 
piles  2  through  5. 

The  final  pile  layout  is  not  the  same  as  in  the  linked 
examples.  Pile  (1)  is  orientated  outward.  The  other  piles 
(2, 3, 4, 5)  are  more  vertical  than  pile  (1).  This  relatively 
small  loading  does  not  require  all  piles  to  be  orientated 
outward.  For  this  magnitude  of  loading  a  better  grouping 
might  be  the  central  piles  in  one  group  and  Pile  (1)  in  second 
group . 
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Table  4-8.  Initial  layouts  and  sizes  for  examples  1.4, 


and  1.5. 


|  Pile 

|  Group 

mm 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

1 

729. 

1. 

180. 

0. 

2 

729. 

5. 

180. 

0. 

3 

729. 

10. 

180. 

0. 

4 

729. 

10. 

0. 

0. 

5 

729. 

1. 

0. 

0. 

Table  4-9.  Final  layouts  and  sizes  for  example  1.4 


Pile 

Group 


(in*) 

Batter 

729. 

1.30 

729. 

3.51 

729. 

3.71 

729. 

4.21 

729. 

5.80 

Phi 

(degrees) 


Theta 

(degrees) 


1 


2 


3 


(a) 


Arrow*  indlcat* 
balt*r  direction. 


2.5’ 


Piles  are  ungrouped. 


depth  ■  30' 
{b} 


Figure  4-10.  Initial  pile  layouts  for  examples  1.4  and  1.5. 
{a}:Plan  view,  {bh-Elevation. 


1 


2 


depth  ■  30’ 
lb} 


Figure  4-11.  Final  pile  layouts  for  example  1.4 
{a}:Plan  view,  {bhElevation. 
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1.5E4 


Weight 

(lbs.) 


1.0E4 


5000.0 


0 


5 


10 


15 


Iteration 


Figure  4- 1  2.  Convergence  of  weight  for  examples  1 .4  and  1 .5. 
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Example  1.5:  Some  of  the  best  weight  reductions  are 
expected  from  the  elimination  of  piles  from  the  design.  The 
elimination  of  a  pile  can  be  done  by  discrete  methods.  A 
continuous  method  for  the  elimination  of  piles  was  attempted. 
It  was  hoped  that  elimination  of  the  minimum  member  size 
constraints  would  allow  a  member  to  reduce  itself  in  size 
until  it  reaches  zero.  In  this  example  there  is  no  minimum 
member  size  constraints.  See  the  initial  layouts  in  Tables  4- 
1,8  and  Figure  4-10. 

See  the  final  layout  in  Table  4-10,  and  Figure  4-13. 
After  15  iterations  the  weight  reduces  to  9360.  lbs  from 
12799.  lbs  as  shown  in  Figure  4-12.  The  weight  increased  in 
the  third  iteration  because  the  first  three  layouts  had  stress 
constraint  violations. 

The  members  have  reduced  in  size  to  a  roughly  equivalent 
value.  It  is  not  obvious  if  a  member  can  be  eliminated.  This 
method  of  member  elimination  is  not  satisfactory  for  this 
example.  It  is  expected  that  for  larger  examples  that  this 
method  will  behave  more  beneficially. 

From  Figure  4-12  it  is  possible  to  achieve  further  weight 
reductions  by  continuing  the  optimization  process.  As  stated 
above  the  method  of  this  example  will  not  be  used  therefore 
continuation  of  the  process  is  unnecessary. 
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The  active  constraints  of  the  final  layout  are:  the  tip 
interferences  for  piles  2  through  5;  the  minimum  batter  for 
pile  1;  and  the  stresses  in  piles  1,5,2,  and  3  for  load  cases 
1,1,2,  and  2  respectively.  Note  that  many  stress  constraints 
are  active.  Since  the  member  sizes  are  not  limited  the 
optimization  algorithm  improves  the  weight  until  other 
constraints  become  active. 


Table  4-10.  Final  layouts  and  sizes  for  example  1.5. 


Pile 

Group 

IHH 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

1 

800.8 

1.00 

180. 

0. 

2 

543.3 

2.34 

0. 

0. 

3 

451.2 

2.42 

0. 

0. 

4 

449.9 

2.63 

0. 

0. 

5 

349.4 

3.18 

0. 

0. 
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depth  a  30’ 
(bj 


Figure  4-13.  Final  pile  layouts  for  example  1.5 
(ahPlan  view,  {b):Elevation. 
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4.2  EXAMPLE  2. 


Example  2.  This  example  has  eight  piles.  None  of  the 
piles  are  linked  to  any  of  the  other  piles.  The  piles  are 
arranged  in  an  unsymmetrical  pattern  as  shown  in  Tables  4- 
11,12  and  Figure  4-14.  The  minimum  size  of  the  cross-sections 
was  not  constrained.  This  example  was  a  second  attempt  to 
eliminate  piles  by  allowing  them  to  reduce  to  zero.  This 
example  has  the  largest  number  of  groups  and  therefore  the 
largest  number  of  variables  of  any  of  the  problems. 

The  final  layout  is  in  Table  4-13  and  Figure  4-15.  As  in 
example  1.5  the  optimization  process  did  not  produce  a  list  of 
piles  which  could  be  obviously  eliminated.  A  few  of  the  piles 
are  slightly  larger  than  the  others.  Piles  (2, 3, 6, 8)  all  have 
similar  moments  of  inertia  and  batter  directions.  Piles 
(1,4, 5, 7)  also  have  similar  properties.  The  piles  may  be 
assigned  to  two  linked  groups. 

The  weight  improved  from  53856.  to  20341.  lbs  after  15 
iterations  as  shown  in  Figure  4-16.  While  developing  the 
optimization  algorithm  it  was  noted  that  a  higher  convergence 
control  parameter  was  required  for  this  example.  This  is 
because  the  cross-section  variables  have  no  restrictions  on 
their  size.  They  are  free  to  greatly  change  in  magnitude  at 
each  iteration.  The  higher  convergence  control  parameter 
restricts  the  variables  from  large  changes  that  can  result  in 
oscillations.  Figure  4-16  shows  that  the  weight  had 
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effectively  converged  after  five  iterations.  After  that  point 
in  any  iteration  the  weight  and  the  variables  do  not  greatly 
change.  The  changes  are  never  small  enough  to  meet  the 
convergence  criteria,  therefore  an  even  larger  convergence 
control  parameter  could  be  used. 

The  active  constraints  of  the  final  layout  are:  the 
stresses  in  piles  l,and  7  for  load  case  1;  and  all  piles  for 
load  case  2.  The  final  layout  is  not  a  feasible  solution. 
Piles  6  and  7  intersect  each  other.  A  tip  interference 
constraint  for  piles  6  and  7  was  not  specified  as  a  constraint 
to  consider. 
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Table  4-11.  Loading  for  example  2. 


Load 

Px 

Py 

Px 

My 

Mx 

case 

kips 

kips 

kips 

ft-k 

ft-k 

ft-k 

1 

-240. 

320. 

-1000. 

0. 

0. 

2 

0. 

0. 

500. 

0. 

0. 

0. 

Table  4-12.  Initial  layouts  and  sizes  for  example  2. 


Pile 

Group 


1, 2,3,4, 
5/ 6,7, 8 


Batter 

Phi 

(degrees) 

Theta 

(degrees) 

729. 

100. 

90. 

90. 

Table  4-13.  Final  layouts  and  sizes  for  example  2. 


Pile 

Group 

■H 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

1 

254.0 

1.82 

270. 

90. 

2 

160.3 

1.36 

90. 

90. 

3 

166.9 

1.39 

90. 

90. 

4 

337.8 

2.17 

270. 

90. 

5 

282.4 

2.87 

270. 

90. 

6 

148.4 

1.15 

90. 

90. 

7 

267.4 

2.01 

270. 

90. 

8 

158.7 

1.26 

90. 

90. 
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Figure  4-14.  Initial  pile  layouts  for  example  2. 

(a):  Plan  view,  (b):  Elevation. 
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5’  I  5’  I  10’  I  5’  15'  I  Piles  are  ungrouped. 


(a) 


(b)  depth  -  90’ 

Figure  4-15.  Final  pile  layouts  for  example  2. 

(a):  Plan  view,  (b):  Elevation. 
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4.3  EXAMPLE  3. 


Example  3 .  This  example  has  29  piles.  The  piles  are 
arranged  in  three  rows.  A  pile  is  linked  to  the  other  piles 
in  the  same  row  as  shown  in  Tables  4-14,15  and  Figure  4-17. 
This  is  a  larger  problem  and  therefore  probably  a  more 
practical  problem  than  those  previously  discussed.  since 
there  are  only  three  groups  the  number  of  variables  is 
relatively  small. 

See  the  final  layout  in  Table  4-16  and  Figure  4-18.  The 
weight  improved  from  237664.  to  154912.  lbs  after  15 
iterations  as  shown  in  Figure  4-19.  The  active  constraints 
were  the  stress  in  pile  11  for  load  case  2;  and  the  minimum 
sizes  for  all  piles.  See  the  convergence  of  the  size  cf  group 
3  in  Figure  4-20.  it  was  found  that  in  this  problem  a  higher 
convergence  control  parameter  was  required  to  achieve 
convergence . 


53 


Table  4-14.  Loading  for  examples  3.  -  3.2. 


Load 

Px 

Py 

Px 

M* 

My 

Mz 

case 

kips 

kips 

kips 

ft-k 

ft-k 

ft-k 

1 

0. 

-1300. 

1500. 

-10000. 

0. 

0. 

2 

0. 

1580. 

0. 

0. 

0. 

Table  4-15.  Initial  layouts  and  sizes  for  example  3.0 


Pile 

Group 

(in'*) 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

1 

1220. 

100. 

90. 

90. 

2 

1220. 

100. 

90. 

90. 

3 

1220. 

100. 

90. 

90. 

Table  4-16.  Final  layouts  and  sizes  for  example  3.0 


I  Pile 

|  Group 

(il?) 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

i 

729. 

3.77 

270. 

90. 

2 

729. 

2.39 

90. 

90. 

3 

729. 

4.16 

270. 

90. 
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cmcmcm3 


Figure  4-17.  Initial  pile  layouts  for  examples  3.  and  3.1. 
(a):  Plan  view,  (b):  Elevation. 
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(b) 


Depth  ■  70’ 


Figure  4-18.  Final  pile  layouts  for  example  3.  and  3.1. 
(a):  Plan  view,  (b):  Elevation. 
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Legend: 

O  —  O  ex  3. 
□  —  □  ex  3.1 
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A; 
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/A 


0 
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Iteration 


10 


15 


Figure  4—19.  Convergence  of  weight  for  examples  3. -3. 2. 


2000 i - 

1500  -- 
Inertia 

(in4)  °-o-o 

1000--  \ 

729  O-M^i^0---O-O-O-O-O-O-O-O-O-O-C) 
500 -■ 

o  4 - 1 - 1 - 

0  5  10  15 

Iteration  number 

Figure  4-20.  Example  3.2,  convergence  of  sire  of  group  3 
to  the  minimum  size. 
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Example  3.1:  This  example  has  a  different  starting  point 
than  example  3  as  shown  in  Tables  4-14,17  and  Figure  4-17. 
The  initial  cross-sections  are  much  smaller.  The  resulting 
initial  stress  constraints  indicate  the  design  has  constraint 
violations.  In  load  case  2  all  of  the  piles  initially  have  a 
41%  stress  constraint  violation. 

See  the  final  layout  in  Table  4-18  and  Figure  4-18. 
After  15  iterations  the  optimal  point  that  is  reached  is 
almost  exactly  the  same  as  in  example  3.  The  weight 
originally  increases  due  to  the  design  starting  in  the 
infeasible  region.  See  Figure  4-19.  Example  3.  and  3.1 
converge  to  the  same  point.  The  optimum  reached  does  not 
appear  to  be  sensitive  to  the  initial  cross-section  selection. 
This  example  also  demonstrates  that  given  a  highly  over¬ 
stressed  initial  design  the  feasible  design  space  can  be 
reached . 
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Table  4-17.  Initial  layouts  and  sizes  for  example  3 . 1 


I  Pile 

|  Group 

mm 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

I  1 

729. 

100. 

90. 

90. 

2 

729. 

100. 

90. 

90. 

3 

729. 

100. 

90. 

90. 

Table  4-18.  Final  layouts  and  sizes  for  example  3.1 


Pile 

Group 

■H 

Batter 

Phi 

(degrees) 

Theta  I 
(degrees)  f 

1 

748.6 

3.96 

270. 

90.  1 

2 

921.0 

2.35 

90. 

90.  1 

3 

856.3 

4.14 

270. 

90.  1 
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Example  3.2:  The  starting  point  of  example  3.  was  varied 
so  that  the  initial  batters  were  opposite  of  their  apparent 
optimal  position  as  shown  in  Tables  4-14,19  and  Figure  4-21. 

See  the  final  layout  in  Table  4-20  and  Figure  4-22.  The 
weight  changes  from  155839.  to  198157.  lbs  in  eight 
iterations.  See  Figure  4-19.  The  weight  increases  due  to  the 
infeasibility  of  the  initial  design.  The  active  constraints 
are  the  minimum  sizes  for  groups  1  and  3;  the  minimum  batters 
for  groups  1  and  3 ;  and  the  stress  in  pile  11  for  load  case  2 . 

This  example  has  a  much  worse  final  design  than  examples 
3.  and  3.1.  This  demonstrates  that  the  design  space  does  have 
local  optimum  points.  Optimization  algorithms  have  not  been 
developed  which  can  guarantee  that  the  global  optimum  of  a 
nonconvex  problem  can  be  found.  A  common  method  to  search  for 
the  optimum  is  to  use  multiple  initial  designs.  Two  possible 
starting  points  exist  for  each  batter.  They  are  the  opposite 
sides.  An  optimization  algorithm  can  be  developed  to  go 
through  each  of  the  possibilities.  At  this  time  it  is 
expected  that  the  experienced  designer  will  select  a  few 
efficient  initial  layouts  then  choose  the  best  final  design 
produced . 
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Table  4-19.  Initial  layouts  and  sizes  for  example  3.2 


Pile 

Group 

mm 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

1  1 

729. 

3. 

90. 

90. 

1  2 

729. 

3. 

270. 

90. 

1  3 

729. 

3. 

90. 

90. 

Table  4-20.  Final  layouts  and  sizes  for  example  3.2 


I  Pile 

1  Group 

wm 

Batter 

Phi 

(degrees) 

.. 

Theta 

(degrees) 

729.0 

1.00 

90. 

90. 

2 

785.2 

2.20 

270. 

90. 

1  r 

729.0 

1.00 

90. 

90. 
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h  t  h  t  h  ti  h  h  h  t 


1 1 


20 


Spacing  tame  as 
example  3.0 


1 

2 

3 


Arrows  indicate 
batter  direction 


10 


19 


(a) 


Groups: 

□  1 


Pile  group  1  Pile  group  2  Pile  group  3 

(b) 

Figure  4-21.  Initial  pile  layouts  for  example  3.2 
(a)  Plan  view,  (b)  Elevation. 


62 


cmim[m 


cT  i 
2 
3 


r pacing  same  as 
example  3.0 


Arrows  indicate 
batter  direction 


(b) 

Figure  4-22.  Final  pile  layouts  for  example  3.2 
(a)  Plan  view,  (b)  Elevation. 
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4.4  EXAMPLE  4. 


Example  4 .  This  example  has  fifteen  piles  arranged  into 
three  linked  groups  as  shown  in  Tables  4-21,22  and  Figure  4- 
23.  Piles  1  through  9  are  linked  into  a  group.  The  other 
groupings  are:  piles  10  through  12,  and  piles  13  through  15. 
Only  the  piles  in  the  first  group  are  aligned  in  the  plane  of 
the  primary  loading.  This  example  is  a  pseudo  three- 
dimensional  optimization  problem  because  the  piles  do  not  vary 
within  the  same  planes.  This  example  is  not  a  fully  three- 
dimensional  because  the  angle  phi  is  not  a  variable. 

See  the  final  layouts  in  Table  4-23  and  Figure  4-24.  The 
weight  converges  from  45120.  to  39878.  lbs  in  ten  iterations 
as  shown  in  Figure  4-25.  The  weight  originally  increases  due 
to  infeasibility  of  the  initial  design.  The  active 
constraints  are  the  minimum  sizes  for  groups  1  and  2;  and  the 
stresses  in  piles  7,10,13,  and  13  for  load  cases  1,1,1  and  2 
respectively. 
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Table  4-21.  Loading  for  example  4. 


Load  Px  Py  Pt 

case  kips  kips  kips 


-250. 


Mx 

ft-k 


My 

ft-k 


200.  300.  650. 


Table  4-22.  Initial  layouts  and  sizes  for  example  4. 


Theta 

(degrees) 


Table  4-23.  Final  layouts  and  sizes  for  example  4. 


Pile 

Group 


Batter 


Phi 

(degrees) 


Theta 

(degrees) 


729.0 

991.3 


1.61 

1.47 


90. 


Figure  4-23.  Initial  pile  layouts  for  example  4. 
(a)  Plan  view,  (b)  Elevation. 
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Spacing  5 


a 


h — i 


Spacing  10*  Y 
(a) 


Arrows 
batter  d 


Groups: 

1  □ 


depth  -  30' 


15 


9  8,7 

(b) 


12 


Figure  4-24.  Final  pile  layouts  for  example  4. 
(a)  Plan  view,  (b)  Elevation. 


indicate 

irection. 
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4.5  EXAMPLE  5. 


Example  5.  This  example  has  240  pile  arranged  in  three 
linked  groups.  This  example  is  taken  from  Example  3  of  the 
CPGA  manual  (9) .  The  initial  layouts  are  found  in  Tables  4- 
24,25  and  Figure  4-26.  The  only  alteration  to  the  example  is 
that  all  pile  tips  reach  a  depth  of  90*.  Note  that  the  pile 
cap  is  sloped. 

See  the  final  layout  in  Table  4-26.  The  weight 
effectively  converges  after  8  iterations  but  the  optimization 
algorithm  is  allowed  to  proceed  for  15  iterations.  The  weight 
improves  from  1477348.  to  1435684.  lbs  as  shown  in  Figure 
4-27.  The  weight  increases  in  the  first  iteration  due  to  the 
linear  representation  of  the  weight  gradients.  The  algorithm 
expected  the  weight  to  decrease.  The  linear  gradients  do  not 
fully  represent  the  design  space.  A  non-benef icial  iteration 
can  be  performed.  If  non-benef icial  iterations  occur  then 
they  usually  occur  early  in  the  optimization  process  while  the 
variables  are  changing  the  greatest.  When  the  variables 
change  greatly  the  actual  weight  has  a  greater  chance  of 
straying  farther  from  the  linearized  prediction. 

The  active  constraints  are  the  minimum  sizes  for  all 
groups;  and  the  stress  in  pile  229. 
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Table  4-24.  Loading  for  example  5. 


Load 

Px 

Py 

Px 

My 

M* 

case 

kips 

kips 

kips 

ft-k 

ft-k 

ft-k 

1 

0. 

-12000. 

20000. 

-1.E6 

0. 

0. 

Table  4-25.  Initial  layouts  and  sizes  for  example  5. 


Pile 

Group 

mm 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

1 

729. 

100. 

90. 

90. 

2 

729. 

2. 

90. 

90. 

3 

729. 

2. 

270. 

90. 

Table  4-26.  Final  layouts  and  sizes  for  example  5. 


Pile 

Group 

(in*) 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

1 

729. 

6.58 

270. 

90. 

2 

729. 

3.96 

90. 

90. 

3 

729. 

3.92 

270. 

90. 
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*2.000* 


5 


10 


Iteration 

:nce  of  Weight  for  example  5. 


4.6  EXAMPLE  6. 


Example  6.  This  example  has  18  piles  arranged  in  two 
groups.  This  example  is  taken  from  Example  5  of  the  CPGA 
manual  (9) .  The  only  alterations  are  that  the  depth  of  piles 
1-12  is  held  at  65' ,  and  the  depth  of  piles  13-18  is  held  at 
74'.  A  displacement  limit  of  0.35  inches  was  arbitrarily 
added.  See  the  initial  layouts  in  Tables  4-27,28  and  Figure 
4-28.  Note  that  the  angle  of  orientation  of  each  of  the  piles 
is  not  the  same  as  the  other  piles  within  the  group.  This  is 
a  pseudo  three-dimensional  example. 

See  the  final  layout  in  Table  4-29.  After  seven 
iterations  the  weight  converges  from  95718.  to  87728.  lbs  as 
shown  in  Figure  4-29.  The  active  constraints  are  the  minimum 
sizes  of  both  groups;  and  the  displacement  component  1.  See 
the  displacement  convergence  in  Figure  4-30. 
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Table  4-27.  Loading  for  example  6. 


Load 

case 

P* 

kips 

V  ,Py 

kips 

P* 

kips 

ft-k 

My 

ft-k 

M* 

ft-k 

1 

-142.8 

0. 

1272. 

0. 

5600. 

0. 

Table  4-28.  Initial  layouts  and  sizes  for  examples  6., 
and  6.1. 


1  Pile 

1  Group 

(in4) 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

i 

904. 

4. 

See  Figure. 

90. 

1  2 

729. 

100. 

See  Figure. 

90. 

Table  4-29.  Final  layouts  aid  t.  zes  for  example  6. 


Group 

(in4) 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

1 

729. 

0 

CO 

• 

CO 

Unchanged . 

90. 

2 

729. 

14.81 

Unchanged . 

90. 
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ELEVATION 


Example  6 . 1  This  example  has  the  same  starting  point  as 
example  6.  The  moment  of  inertias  are  held  constant  to 
demonstrate  the  final  phases  of  design.  Discrete  section 
sizes  can  be  selected  and  held  fixed.  The  example  is  then 
optimized  with  respect  to  the  batters. 

See  the  final  layout  in  Table  4-30.  The  weight  improves 
to  94830.  after  seven  iterations  as  shown  in  Figure  4-29. 
Holding  the  inertias  constant  causes  the  final  weight  to  be 
much  higher  than  when  they  are  allowed  to  vary.  The  only 
active  constraint  is  the  displacement  component  1.  See  the 
convergence  of  the  displacement  component  in  Figure  4-31. 

Table  4-30.  Final  layouts  and  sizes  for  example  6.1. 


Pile 

Group 

(in4) 

Batter 

Phi 

(degrees) 

Theta 

(degrees) 

1 

904. 

9.13 

Unchanged . 

90. 

2 

729. 

17.27 

Unchanged . 

90. 

Table  4-31.  Example  soil  type  and  fixity. 


Example 


Soil  type 

Pile  Fixity 

NH  0.019 

All  fixed 

NH  0.019 

All  pinned 

NH  0.019 

All  pinned 

NH  0.019 

All  fixed 

NH  0.019 

All  pinned 

ES  0.312 

All  pinned 
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0  1  2  3  4  5  6  7 

Iteration 

Figure  4-31.  Satisfaction  of  displacement  limit  for  example  6.1 
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5.  SUMMARY  AND  CONCLUSIONS 

An  optimality  criteria  approach  was  adapted  to 
structurally  optimize  the  weight  of  steel  HP-14  piles  used  to 
support  structures  which  are  modeled  as  rigid  slabs.  The 
examples  included  dams,  a  retaining  wall,  and  a  control  tower. 
The  optimization  method  is  an  iterative  numerical  process.  An 
optimality  approach  was  implemented  because  it  converges 
quickly  to  an  optimal  solution. 

The  primary  strategy  of  the  optimality  criteria  approach 
is  to  form  an  optimality  criteria  by  using  the  Lagrange 
equation  for  the  weight  and  constraints.  The  optimality 
criteria  is  used  to  form  a  variable  recurrence  equation. 

The  secondary  problem  is  to  solve  for  the  Lagrange 
Multipliers  by  anticipating  which  constraints  are  active  and 
then  solving  a  set  of  linear  equations.  The  algorithm 
presented  for  selection  of  active  constraints  has  worked  for 
all  examples,  and  for  all  other  feasible  problems  that  were 
attempted.  The  selection  process  uses  a  large  portion  of  the 
computation  time  therefore  this  is  an  area  of  ootential 
research. 

A  method  was  developed  to  promote  the  best  convergence  of 

topological  variables  with  optimality  criteria.  The 

coordinate  system  for  each  topological  variable  is  selected 

according  to  the  predictability  of  the  variable,  and  according 

to  how  the  weight  changes  with  the  variable.  Since 
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topological  variables  can  have  a  negative  weight  gradient, 
the  optimality  criteria  methods  were  adapted  to  properly 
handle  such  variables  by  altering  the  coordinate  systems. 

An  alternative  energy  objective  function  is  used  for  a 
variable  when  the  variable  has  no  direct  effect  on  the  weight. 
A  weight  gradient  is  required  for  all  variables  with  the 
current  optimality  criteria  version.  Energy  is  chosen  to  be 
minimized  because  it  promotes  lower  displacements  which  are 
generally  desirable. 

The  pile  cross-sections  and  pile  batters  were  allowed  to 
vary.  The  variables  were  not  held  to  discrete  values  but  were 
allowed  to  change  in  a  continuous  manner.  The  pile  cross- 
sections  were  represented  by  a  single  primary  variable  which 
was  the  section's  major  axis  inertia.  Equations  were 
developed  to  closely  approximate  the  other  cross-section 
variables  such  as  area  given  the  inertia. 

Six  primary  examples  were  optimized.  The  number  of  piles 
in  each  of  the  examples  ranged  from  5  to  240  piles.  The  piles 
were  either  placed  in  linked  groups  or  they  were  allowed  to 
independently  vary.  Linking  piles  into  groups  provides 
consistency  that  is  desired  when  installing  the  piles,  and 
linking  eases  computational  requirements. 

The  designs  which  had  initial  constraint  violations 

quickly  approached  the  feasible  design  space.  The  initially 

feasible  problems  had  slight  to  significant  weight  improvement 
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from  the  initial  design.  The  number  of  iterations  for  the 
designs  to  converge  was  relatively  small.  In  a  few  iterations 
the  design  was  usually  within  the  neighborhood  of  the  final 
converged  design. 

As  a  few  of  the  designs  approached  the  converged  design 
the  variables  and  weight  oscillated.  The  estimate  of  the 
optimum  point  using  the  linearizations  of  the  constraints  and 
weight  function  can  overshoot  or  undershoot  the  optimum. 
Overshooting  and  the  resulting  oscillation  can  be  reduced  by 
increasing  the  convergence  control  parameter.  The 
linearizations  also  caused  Example  5  to  have  a  non-benef icial 
iteration.  The  linear  algorithm  expected  a  weight  improvement 
in  the  first  iteration  even  though  it  did  not  occur. 

Locally  optimal  points  with  weights  greater  than  the 
global  optimum  may  be  reached.  It  was  demonstrated  that  the 
optimization  algorithm  can  be  sensitive  to  the  choice  of  the 
initial  batters,  but  it  does  not  appear  as  sensitive  to  the 
choice  of  initial  cross-sections.  The  experienced  designer 
should  use  multiple  starting  points  for  the  optimization 
process.  The  choice  of  which  piles  to  link  will  have  an 
effect  on  the  final  optimal  solution.  The  designer  may  choose 
to  leave  the  piles  unlinked  until  the  optimum  is  approached. 
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7.  APPENDIX:  PILE  GROUP  BEHAVIOR 


The  optimization  methods  are  applied  to  an  existing  pile 
group  analysis  computer  program.  The  analysis  methods  are  not 
altered  by  the  optimization  algorithm.  The  following  is  a 
discussion  of  the  pile  group  analysis  methods  which  is  taken 
from  BASIC  PILE  GROUP  BEHAVIOR  (10),  pp.  6,7. 

The  "basic  method  of  pile  group  analysis  ...  is  valid 
for  static  analysis  of  a  linear,  elastic  system.  Interaction 
between  pile  and  structure  is  limited  to  the  extremes  of  a 
fully  fixed  or  fully  pinned  connection.  Interaction  between 
the  pile  and  soil  is  represented  by  a  linear,  elastic  pile 
stiffness  (applied  load  per  unit  deflection)  at  the  top  of  the 
pile.  The  base  of  the  structure  is  assumed  to  act  as  a  rigid 
body  pile  cap  connecting  all  piles;  the  cap  flexibility  is  not 
considered.  This  method  of  analysis  will  also  be  incorporated 
in  the  new  pile  group  analysis  program. 

"The  basic  pile  group  analysis  method  represents  each 
pile  by  its  calculated  stiffness  coefficient,  in  the  manner 
proposed  by  Saul.  The  stiffness  coefficients  of  all  piles  are 
summed  to  determine  a  stiffness  matrix  for  the  total  pile 
group.  Displacements  of  the  rigid  pile  cap  are  determined  by 
multiplying  the  sets  of  applied  loads  by  the  inverse  of  the 
group  stiffness  coefficients  to  determine  the  forces  acting  on 
each  pile  head.  The  key  step  in  the  method  is  in  determining 
individual  pile  stiffness  coefficients,  at  the  pile  head, 
based  on  known  or  assumed  properties  of  pile  and  soil.  Since 
this  is  a  three-dimensional  analysis  method,  each  pile  head 
has  six  degrees  of  freedom  (DOF) ,  three  translations  and  three 
rotations.  A  stiffness  coefficient  must  be  determined  for 
each  DOF  and  for  all  coupling  effects  (e.g.  lateral  deflection 
due  to  applied  moment) .  The  pile  location  and  batter  angle 
are  also  accounted  for  when  individual  pile  stiffness 
coefficients  are  combined  to  form  the  total  stiffness  matrix 
for  the  pile  group. 

"Piles  are  mathematically  represented  in  the  analysis  by 
their  axial,  lateral  and  rotational  stiffness,  as  springs 
resisting  motion  of  the  rigid  cap." 

The  pile  stresses  are  evaluated  to  determine  the 
feasibility  of  the  design  by  use  of  interaction  equations. 
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The  following  equations  are  taken  from  BASIC  PILE  GROUP 
BEHAVIOR  (10),  pp.  A8,A9. 

"Steel  piles  subject  to  axial  load  and  bending  shall  be 
proportioned  to  satisfy  the  following  requirements: 

£a  +  Cmx^by  +  ^my^by  ^  ^  _  q 

Fa  (l-fa/F'ax)  (1  -fjF'ey)Fby 


and: 

—  +  i  1,0  ( when —  s  0.15) 

Fa  FbX  Fby  Fa 


where : 

pt  = _ k2E _ 

8  F.S.  {KtLjxJ2 


and: 


fa 


computed  axial  stress  (psi) 


fta  or  fby  =  computed  compressive  bending  stress  about 
the  x  axis  and  y  axis,  respectively  (psi) 


Fa  =  allowable  axial  stress  (psi) 

Fbx  or  Fby  =  allowable  compressive  bending  stess  about 
the  x  and  y  axis,  respectively  (psi) 


E  = 
Lb  = 

Kb  - 

rb  = 


C«*  or 


my 


F.S.  = 


modulus  of  elasticity  (29,000,000  psi) 
actual  unbraced  length  of  pile  in  the 
plane  of  bending  (inches) 
effective  length  factor  as  defined  by 
AISC  in  the  plane  of  bending  (inches) 
radius  of  gyration  in  the  plane  of 
bending  (inches) 

coefficient  about  x  and  y  axes, 
respectively,  as  defined  by  AISC 
Factor  of  Safety 


The  lower  regions  of  the  bearing  piles  are  subject  to 


axial  stresses.  Since  the  lower  region  of  the  pile  is  subject 
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to  damage  during  driving  a  greater  factor  of  safety  is  used. 

The  axial  load  at  the  tip  must  be  less  than  the  damaged  cross- 

section  capacity  and  the  soil  bearing  capacity. 

The  point  of  maximum  moment  in  a  pile  is  dependant  upon 

the  soil  conditions  and  the  fixity  of  the  pile  to  the  cap. 

The  following  is  taken  from  the  CPGA  Manual  (9),  pp.  44,45. 

"The  pile  forces,  used  for  the  above  comparisons  with  the 
allowables,  are  those  forces  calculated  as  acting  on  the  top 
of  the  pile.  "Pinned"  piles,  however,  have  no  moment  at  the 
top,  but  do  experience  bending  beneath  the  top,  depending  on 
lateral  loads  and  soil  properties.  The  critical  moments  may 
be  approximated  for  design  purposes  as  a  constant  times  the 
lateral  force  on  the  pile 

Ml  =  KMPl  F2 

where 

Ml  =  design  moment 
F2  =  lateral  force 
KMPl  =  constant 

"These  calculated  moments  are  then  used  in  the  allowable  load 
comparisons.  When  pile  and  soil  properties  have  been  given 
for  a  pinned  pile,  the  design  moment  factors,  KMPl  and  KMP2 , 
may  be  calculate  automatically." 
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8 .  NOMENCLATURE 


WT  =  total  weight  of  the  steel  piles. 

^  =  density  of  steel. 

V;  =  volume  of  element  i. 

N  =  number  of  piles, 

hj  =  constraint  j  . 

m  =  number  of  constraints. 

Oj  =  stress  in  member  j . 

=  upper  limit  on  the  stress  in  member  j . 

I„  =  major  axis  moment  of  inertia. 

Ixx  -  lower  bound  on  the  major  axis  moment  of  inertia 
I]C-  =  upper  bound  on  the  major  axis  moment  of  inertia 
L  =  Lagrangian. 

\j  =  Lagrange  Multiplier  for  the  j*  constraint, 

dj  =  design  variable  number  i. 

n  =  number  of  variables, 

r  =  convergence  control  parameter, 

n,  =  number  of  active  variables. 

Iyy  =  minor  axis  moment  of  inertia  of  a  given  member. 
A  =  area  of  a  given  member. 

c,  =  extreme  fiber  distance  in  the  direction  of  a 
member's  local  x  axis. 
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cy  =  extreme  fiber  distance  in  the  y  direction. 
Pi  =  global  force  component  i. 

D;  *  global  displacement  component  i. 
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straint  Processing,  Volumes  I  and  II 

A  Case  Committee  Study  of  Finite  Element  Analysis  of  Concrete 
Flat  Slabs 

User's  Guide:  Computer  Program  for  Two-Dimensional  Analysis 
of  U-Frame  Structures  (CUFRAM) 

User's  Guide:  For  Concrete  Strength  Investigation  and  Design 
(CASTR)  in  Accordance  with  ACI  318-83 

Finite-Element  Method  Package  for  Solving  Steady-State  Seepage 
Problems 

User's  Guide:  A  Three  Dimensional  Stability  Analysis/Design 
Program  (3DSAD)  Module 
Report  1 :  Revision  1 :  General  Geometry 
Report  2:  General  Loads  Module 
Report  8:  Free-Body  Module 
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Date 


Instruction  Report  ITL-87-4 
Technical  Report  ITL-87-4 


Instruction  Report  GL-87-1 

Instruction  Report  ITL-87-5 
Instruction  Report  ITL-87-6 

Technical  Report  ITL-87-8 

Instruction  Report  ITL-88-1 

Technical  Report  ITL-88-1 

Technical  Report  ITL-88-2 

Instruction  Report  ITL-88-2 

Instruction  Report  ITL-88-4 

Instruction  Report  GL-87-1 

Technical  Report  ITL-89-3 
Technical  Report  ITL-89-4 


User's  Guide:  2-D  Frame  Analysis  Link  Program  (LINK2D) 

Finite  Element  Studies  of  a  Horizontally  Framed  Miter  Gate 

Report  1 :  Initial  and  Refined  Finite  Element  Models  (Phases 
A,  B,  and  C),  Volumes  I  and  II 
Report  2:  Simplified  Frame  Model  (Phase  D) 

Report  3:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies-Open  Section 

Report  4:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies-Closed  Sections 

Report  5:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies-Additional  Closed  Sections 
Report  6:  Elastic  Buckling  of  Girders  in  Horizontally  Framed 
Miter  Gates 

Report  7:  Application  and  Summary 

User's  Guide:  UTEXAS2  Slope-Stability  Package:  Volume  I, 
User's  Manual 


Sliding  Stability  of  Concrete  Structures  (CSLIDE) 

Criteria  Specifications  for  and  Validation  of  a  Computer  Program 
for  the  Design  or  Investigation  of  Horizontally  Framed  Miter 
Gates  (CMITER) 

Procedure  for  Static  Analysis  of  Gravity  Dams  Using  the  Finite 
Element  Method  -  Phase  la 

User's  Guide:  Computer  Program  for  Analysis  of  Planar  Grid 
Structures  (CGRID) 

Development  of  Design  Formulas  for  Ribbed  Mat  Foundations 
on  Expansive  Soils 

User's  Guide:  Pile  Group  Graphics  Display  (CPGG)  Post¬ 
processor  to  CPGA  Program 

User's  Guide  for  Design  and  Investigation  of  Horizontally  Framed 
Miter  Gates  (CMITER) 

User's  Guide  for  Revised  Computer  Program  to  Calculate  Shear, 
Moment,  and  Thrust  (CSMT) 

User’s  Guide:  UTEXAS2  Slope-Stability  Package;  Volume  II, 
Theory 

User's  Guide:  Pile  Group  Analysis  (CPGA)  Computer  Group 

CBASIN-Structural  Design  of  Saint  Anthony  Falls  Stilling  Basins 
According  to  Corps  of  Engineers  Criteria  for  Hydraulic 
Structures:  Computer  Program  X0098 
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Technical  Report  ITL-89-5 

Technical  Report  ITL-89-6 
Contract  Report  ITL-89-1 
Instruction  Report  ITL-90-1 
Technical  Report  ITL-90-3 

Instruction  Report  ITL-90-6 
Instruction  Report  ITL-90-2 
Technical  Report  ITL-91-3 

Instruction  Report  ITL-91-1 

Instruction  Report  ITL-87-2 
(Revised) 

Technical  Report  ITL-92-2 
Technical  Report  ITL-92-4 
Instruction  Report  ITL-92-3 

Instruction  Report  ITL-92-4 
Instruction  Report  ITL-92-5 
Contract  Report  ITL-92-1 
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Title 

CCHAN-Structural  Design  of  Rectangular  Channels  According 
to  Corps  of  Engineers  Criteria  for  Hydraulic 
Structures;  Computer  Program  X0097 

The  Response-Spectrum  Dynamic  Analysis  of  Gravity  Dams  Using 
the  Finite  Element  Method;  Phase  II 

State  of  the  Art  on  Expert  Systems  Applications  in  Design, 
Construction,  and  Maintenance  of  Structures 

User's  Guide:  Computer  Program  for  Design  and  Analysis 
of  Sheet  Pile  Walls  by  Classical  Methods  (CWALSHT) 

Investigation  and  Design  of  U-Frame  Structures  Using 
Program  CUFRBC 

Volume  A:  Program  Criteria  and  Documentation 
Volume  B;  User's  Guide  for  Basins 
Volume  C:  User’s  Guide  for  Channels 

User's  Guide:  Computer  Program  for  Two-Dimensional  Analysis 
of  U-Frame  or  W-Frame  Structures  (CWFRAM) 

User’s  Guide:  Pile  Group-Concrete  Pile  Analysis  Program 
(CPGC)  Preprocessor  to  CPGA  Program 

Application  of  Finite  Element,  Grid  Generation,  and  Scientific 
Visualization  Techniques  to  2-D  and  3-D  Seepage  and 
Groundwater  Modeling 

User's  Guide:  Computer  Program  for  Design  and  Analysis 
of  Sheet-Pile  Walls  by  Classical  Methods  (CWALSHT) 

Including  Rowe's  Moment  Reduction 

User's  Guide  for  Concrete  Strength  Investigation  and  Design 
(CASTR)  in  Accordance  with  ACI 318-89 

Fiinite  Element  Modeling  of  Welded  Thick  Plates  for  Bonneville 
Navigation  Lock 

Introduction  to  the  Computation  of  Response  Spectrum  for 
Earthquake  Loading 

Concept  Design  Example,  Computer  Aided  Structural 

Modeling  (CASM) 

Report  1 :  Scheme  A 
Report  2:  Scheme  B 
Report  3:  Scheme  C 

User's  Guide:  Computer-Aided  Structural  Modeling 

(CASM)  -  Version  3.00 

Tutorial  Guide:  Computer-Aided  Structural  Modeling 

(CASM)  -  Version  3.00 

Optimization  of  Steel  Pile  Foundtions  Using  Optimality  Criteria 
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